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Abstract 

Arecibo S-band (A = 13 cm) radar observations of Comet C/2001 A2 (LINEAR) on 
2001 July 7-9 showed a strong echo from large coma grains. This echo was significantly 
depolarized. This is the first firm detection of depolarization in a grain-coma radar echo and 
indicates that the largest grains are at least A/2 or 2 cm in radius. The grains are moving at 
tens of ms 4 with respect to the nucleus. The non-detection of the nucleus places an upper 
limit of 3 km on its diameter. The broad, asymmetric echo power spectrum suggests a fan 
of grains that have a steep (differential number ~ a -4 ) size distribution at cm-scales, though 
the observed fragmentation of this comet complicates that picture. 

Keywords: radar, comets 
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1. Introduction 

When comets enter the inner solar system, we get a close look at material that is 
relatively pristine, closely resembling the volatile-rich material that formed the outer planets 
and their satellites. Comets may also have contributed significant volatiles to the early 
terrestrial planets during accretion, although the relative contribution of cometary impactors 
and more rocky asteroidal bodies is still debated. Comets also present an impact hazard for 
the Earth. Understanding the physical structure of cometary nuclei and how they fragment 
is important for learning about their formation and evaluating the impact hazard. The 
break-up of Comet D/1993 F2 (Shoemaker-Levy 9) has reinforced the idea that comets may 
be aggregates of icy grains that easily fragment into similar-sized pieces (e.g., Asphaug and 
Benz, 1996). 

The existence of large icy grains in cometary comae has long been inferred. Solar 
radiation pressure affects the orbital motion of small particles (Burns et al., 1979). By 
examining the orbital evolution of grains, it is possible to determine their size, based on 
the extent to which solar radiation affects their orbits. The observed motions of particles in 
comet tails are consistent with sizes of mm to cm in a number of cases. Modeling of comet 
tails suggests large particles in some cases (summarized in Fulle, 2004). The orbital motion 
of some of dust bands also require that the particles be mm to cm in size (e.g., Sykes et al., 
1986, Reach et al., 2000). 

Submillimeter observations of C/1996 B2 (Hyakutake) showed evidence of emission from 
mm-scale particles (Jewitt and Matthews, 1997). The distributions of CO and formaldehyde 
are often found to be inconsistent with a purely nuclear source, and sublimation of grains 
is inferred (e.g., Eberhardt, 1999, Gunnarson, 2003). However, it is not known whether 
these grains are an important source of other gas species or dust. Radar observations have 
indicated the presence of large grains in comet Halley (Campbell et al., 1983) and Iras- 
Araki-Alcock (Harmon et al., 1983), though the grains were not directly detected. Radar 
observation of large (cm-sized) grains near the nucleus determines their velocity distribution, 
improving our understanding of how grains leave the nucleus, which in turn helps determine 
the relative contributions of the nucleus and the grains to the gas coma. 

This paper will give a brief description of earlier radar results (Section 2), present our 
observations of C/2001 A2, and use those observations to constrain the size of the nucleus and 
the population of coma grains (Section 3). We then interpret the grain population both in 
the context of a simple gas-drag ejection model, yielding a size distribution of the coma grains 
that contribute to the radar echo (Section 4), and as the outcome of earlier fragmentation 
of the nucleus (Section 5). We compare the results to those from other comets and using 
other methods (Section 6). We then conclude (Section 7) that the grain size distribution 
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determined in the modeling is likely to be correct even if the model is oversimplified. 



2. Background 

If a comet nucleus comes close enough to the Earth, we may have the opportunity to 
obtain 2-dimensional radar delay-Doppler images, which can provide a wealth of information 
about the comet's shape, size, and surface characteristics. However, because the signal-to- 
noise ratio (SNR) of radar observations is very sensitive to the distance A to the target 
(SNR ~ A~ 4 ), radar imaging of comets requires an approach distance of 0.1 AU or less. Over 
the last twenty years, only comets C/1983 HI (IRAS-Araki-Alcock) and C/1983 Jl (Sugano- 
Saigusa- Fujikawa) in 1983 and comet C/1996 B2 (Hyakutake) in 1996 have approached 
within this distance from the Earth, but radar imaging was not technically feasible for any 
of those objects. 

Even if imaging is not possible, we can measure the basic cm-wavelength scattering 
properties and estimate the size of the nucleus of any comet that comes within about 0.5 
AU of the Earth. Altogether, nine comets have been observed with radar, all at Arecibo 
except for IRAS-Araki-Alcock, which was observed at both Arecibo and Goldstone, and 
Hyakutake, which was observed only at Goldstone. All of these objects were observed by 
transmitting a monochromatic signal and measuring the Doppler shift and broadening of the 
return echoes. This broadening can yield rotational information for the nucleus and particle 
velocities for coma grains. We measure both circular polarizations of the returned echo, 
because the polarization properties of the echo can yield surface roughness information for 
the nucleus and information on coma grain sizes. Of the nine comets successfully observed, 
the nucleus was detected in six (all but Halley, C/2001 A2, and C/2002 06) and coma grains 
were detected in five (Harmon et al., 2004). 

For a cometary echo, two different mechanisms introduce Doppler broadening. These 
are illustrated in Fig. 1, which shows the echo power spectrum of IRAS-Araki-Alcock taken 
at Arecibo in 1983 (Harmon et al., 1989). The central spike is an echo from the solid comet 
nucleus. Its frequency width 8v is proportional to the rotation rate and diameter of the 
nucleus: 

. Sv AnDusini AnDsini 

ov = 2 — v = = , 1 

c Pc PA ' K J 

where 5v is the velocity difference, with respect to the observer, between the approaching and 
receding edges, P is the rotation period, D is the nucleus equatorial diameter in the plane of 
the sky, i is the angle between the comet's spin axis and the line of sight, v is the observing 
frequency, A is the observing wavelength, and c is the speed of light. The broader echo 
component from grains in the coma and its frequency width is due to the velocity dispersion 
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(projected along the line of sight) of all coma particles within the beam that contribute to 
this part of the radar echo. 

For IRAS-Araki-Alcock, the reflected power from the coma was almost entirely (98.6%) 
in the opposite circular polarization (OC) from that transmitted (see Fig. 1), suggesting 
single scattering from small (Rayleigh) scatterers no larger than about 2 cm in size, but 
not much smaller, either (Harmon et al., 1989). The echo from the nucleus still shows most 
(90%) of the echo in the OC, with the other 10% in the same sense of circular polarization 
(SC) as that transmitted. Multiple scattering or reflection from a surface that is rough at 
wavelength scales has "depolarized" a portion of the echo power. The IRAS-Araki-Alcock 
observations were made when the comet was very close (0.033 AU) to the Earth, resulting 
in extremely strong echoes. No comparable opportunity has presented itself since. 

3. Observations of C/2001 A2 (LINEAR) 

Comet C/2001 A2 (LINEAR) was discovered on 2001 January 3 by the LINEAR near- 
Earth asteroid discovery program (Stokes et al., 2000). We observed the comet on 2001 
July 7-10 with the Arecibo S-band radar. The observing circumstances are summarized in 
Table 1. 

We transmitted a left-circularly-polarized, monochromatic 2380-MHz radio wave at a 
power of 900 kW, and observed the echo power spectrum in both left (SC) and right (OC) 
circular polarizations. We report the radar cross section a, defined as the cross sectional 
area nR 2 of a metal sphere that would reflect the same power as we measure. The radar 
albedo is then the ratio of the radar cross section to the physical cross section of the target. 
The polarization ratio /i c is the ratio o- sc /a oc . 

At the time of observations, the Arecibo telescope had a forward gain of 73 dB. The 
receiving system had an effective collecting area of 9 K/Jy (25000 m 2 ) and a system noise 
temperature of about 25K. The gain, effective area, and system temperature vary by about 
20% with telescope pointing, and we used values based on the actual pointing in data re- 
duction. Every 10 s, we changed the transmitter frequency by 10 kHz, with a corresponding 
change applied in the data processing, in order to obtain a measurement of the background 
so that it could be accurately subtracted. 

It is not feasible to simultaneously transmit and receive with the same antenna, so we 
transmitted for the round-trip light time to the comet and then, shortly before the reflected 
light arrived back at the antenna, turned off the transmitter and received the returned echo 
for the round-trip light time to the comet. The radial component of the overall (ephemeris) 
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velocity of the comet was removed in the datataking process by introducing small offsets in 
the frequency of the transmitted signal to compensate for the Doppler shift of the center of 
mass of the target as predicted by an ephemeris based on previous optical observations. 

C/2001 A2 fragmented several times before (and perhaps after) our observations. For 
these observations, the pointing and Doppler ephemerides were for the "B" fragment, at that 
time apparently the only visible fragment (Sekanina et al., 2002), though other fragments 
and their remnants may have been in the beam, particularly fragments D, E, and F (see 
Section 5). 

Echo power spectra of C/2001 A2 on 2001 July 7-9 are shown in Fig. 2, and their 
parameters are summarized in Table 1. The relatively broad OC echo power spectrum 
peaked near the predicted ephemeris Doppler shift, and had a width of at least 900 Hz (57 
ms" 1 ) on each date. The echo is asymmetric, with more power in the receding (negative 
Doppler shift) direction. Asymmetry was also clear in the coma echoes from IRAS-Araki- 
Alcock (Fig. 1) and Hyakutake (Harmon et al., 1997). 

The SC echo power is about 25% that of the OC power. Its spectrum appears to have a 
somewhat narrower bandwidth than that of the OC echo, although this apparent difference 
may be due to the lower SNR. 

3.1. Measurement Uncertainty 

The uncertainties in the reported values are dominated by calibration errors, not random 
noise: The thermal noise is < 5% for a oc on each day. The systematic errors in the cross 
sections include calibration uncertainties for the transmitter power, the antenna performance, 
and the receiver thermal noise, and total to about 30%. In the polarization ratio /x c , all but 
the receiver noise calibration uncertainties cancel, and the systematic uncertainty is only 
about 5%. The day-to-day variations in these systematic calibration errors are also smaller, 
about 5-10%. The day-to-day variation in the OC and total cross section is 20 — 30%, and 
thus statistically significant, but the other parameters are not. 

3.2. Nucleus or Coma? 

Our conclusion depends on the observed SC echo being reflection from coma grains, so 
we first examine the hypothesis that the echo could come from the nucleus instead. These 
spectra do not show the "spike plus skirt" pattern that distinguished the coma from the 
nucleus echo in (for example) the IRAS-Araki-Alcock spectra. The SC echo appears to be 
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at least 300 Hz wide. The rotation period of the nucleus has been estimated at either 6 or 
3 hours by Woodney et al. (2001). Assuming the same radar albedo as IRAS-Araki-Alcock 
(4%, the lowest well-determined comet radar albedo (Harmon et al., 1999)), a nucleus larger 
than 3 km diameter with rotation period of at least three hours would have been detectable 
with a few Hz bandwidth at > 10a on July 7, The strongest narrow features that could be 
from the nucleus are present only at the 3-a significance level, so we consider 3 km to be an 
upper limit to the nucleus diameter. In order to have a bandwidth of 300 Hz, the nucleus 
would have to be at least 30 km in diameter, which would have been visible optically about 5 
magnitudes brighter than was actually observed in January and February 2001. Thus, from 
the large bandwidth and overall echo spectrum shape we conclude that the echo power was 
dominated by reflections from coma grains. The nucleus may contribute to the peak of the 
echo, but it cannot be distinguished. 



3.3. Coma Echo 



The spectra in Fig. 2 show that there is a substantial SC echo. Although this degree 
of depolarization is typical of comet nuclei, we have shown that the SC echoes cannot have 
come from the nucleus. Hence, we conclude that significant depolarization (/i c = 28%) was 
introduced in the radar backscatter from large-grain coma. This is in stark contrast to 
the coma echo from IRAS-Araki-Alcock (Fig. 1), which had yU c only 1.4%. Cometary grain 
comae are too sparsely populated to give significant depolarization from multiple scattering 
(Harmon et al., 1989, 2004), so any detectable depolarization must arise from the single 
scattering off individual irregular grains. Depolarization from irregular grains is small (~ 
0.1-3%) for sizes in the Rayleigh (a < \/2ir) regime and increases rapidly to large values 
(~ 10-50%) as the radius increases beyond the \/2ir threshold (Harmon et al., 1989). The 
precise depolarization behavior depends on the grain density and degree of irregularity. For 
example, the transition size can be somewhat larger than X/2ir for low-density grains. If the 
differential grain size distribution is described by a power law N a oc a~ a with a cutoff radius 
a max , then the strong Rayleigh dependence for the grain backscatter efficiency will tend to 
make the echo depolarization be dominated by the larger grains unless the size distribution is 
incredibly steep (a > 5), so that ji c can be quite sensitive to whether a max is larger or smaller 
than \/2n. For IRAS-Araki-Alcock, the case was made, based on total-mass arguments and 
the low depolarization, that a max must have been of the order of, but not much smaller than, 
A/27T, which is 2 cm for the 13 cm wavelength used (Harmon et al., 1989), and that a cutoff 
size of this order was consistent with the gravity-limit against gas-drag (see Section 4). A 
comet with a smaller nucleus size or higher gas flux than IRAS-Araki-Alcock could have 
a max larger than \/2tt, and hence show significant depolarization. Both Halley (Campbell 
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et al., 1989) and Hyakutake (Harmon et al., 1997) did, in fact, show hints of non-negligible 
depolarization, although in both cases the SNR was too low to be certain. Thus, the new 
depolarization results for C/2001 A2 provide the first unambiguous radar evidence for large 
grains with sizes of several centimeters or more. 

If the SC spectrum really is narrower than the OC spectrum, it may be due to a 
dispersion in ejection velocity as a function of grain size, because only large particles (> 2 
cm) give substantial SC echo, and larger particles are likely slower moving (see Section 4). 
Note that, unlike the case of a rotating solid body, there is no direct spatial information in 
the Doppler spectrum of the coma: The coma is assumed to fill the radar beam, and we 
observe the echo of the entire population of grains. 

The radar cross-section shows day-to-day variations at the 20% level. These variations 
probably represent real differences in the grain population. There is also a hint of a "shoul- 
der" on the low-frequency tail of the spectrum. Since small grains are more easily accelerated 
to higher velocities than large ones, it is possible that this shoulder appears where the grain 
size approaches the Rayleigh scattering limit (2 cm for these observations), reducing the 
scattering efficiency, and thus the cross section at high velocity. The spectrum of C/2002 06 
(SWAN) shows a similar shoulder (Harmon et al., 2004), but in both cases, the data are too 
noisy to address this idea in detail. 

4. Grain-Coma Modeling 

4.1. Grain Ejection 

We have modeled the grain coma in order to estimate the characteristic grain velocity, 
place constraints on the grain ejection direction, and estimate the nucleus mass-loss rate in 
grains. We adopt the approach taken in modeling the coma radar echoes from comets IRAS- 
Araki-Alcock (Harmon et al., 1989) and Hyakutake (Harmon et al., 1997). The starting 
point is the canonical gas-drag model for dust release originally proposed by Whipple (1951) 
and adapted later by others (see Harmon et al., 1989, Harmon et al., 2004, and references 
therein). 

According to this theory, the radial gas drag from an outgassing nucleus of radius R 
accelerates a grain of radius a from the surface to a terminal velocity V t given by 



V t (a)=C v a-^ 2 (l-a/a m f 2 



(2) 
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where 



C v = 



( 



3C D V g ZR\ 
4p g ) 



(3) 



is a velocity scale factor and 



a 
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is the radius of the largest grain that can be lifted off the surface against nucleus gravity. 
Here, C D is a drag coefficient, V g — 1.7 x 10 3 T(Kelvin) 1 / 2 cms" 1 is the gas velocity, assumed 
to be the thermal expansion velocity at temperature T, Z is the gas mass flux at the surface, 
G is the gravitational constant, p n is the mean density of the nucleus, and p g is the density 
of a grain. 

Since grains are presumably accelerated in jets, where the gas velocity is considerably 
higher than the simple thermal expansion velocity, the scale factor C v will tend to be larger 
for comets with more violent jetting. 

We assume simple forms for the grain ejection directions and grain size distributions. 
The grains are ejected in a conical fan, and the grain radii are assumed to be power-law 
distributed as a~ a . Again, the grains are presumably ejected from illuminated jets, not a 
uniform cone. We compute the Doppler spectrum from the radial components of the terminal 
velocities of the entire ensemble of grains within the radar beam, assuming a stable fan of 
continuously ejected grains, and accounting for their orbital motion. We then adjust C v and 
fan direction to match the width and asymmetry of the observed spectrum. 

In broad terms, the size distribution of the particles determines the shape of the spec- 
trum: a steep distribution will have many small grains, which can be accelerated to a high 
velocity, and give a broad "skirt" to the spectrum. The velocity scale factor C v determines 
the actual width of the spectrum, as long as the size distribution is steep enough that the 
small particles dominate the echo. A fan (as opposed to isotropic) shape is required to give 
an asymmetric spectrum. Changing the cone angle of the fan affects the spectrum width, 
with a narrower cone angle tending to reduce the breadth of the spectrum. However, we find 
that we need a broad cone angle (120°), a steep size distribution and a large C v to model 
the rather broad observed spectrum, so that the effects of adjusting the cone angle will be 
minor, and cannot be easily distinguished from those of adjusting the jet orientation. 

The numerical results that follow assume that the grain density p g = 0.5 gem" 3 and the 
nucleus density p n = lgcm" 3 , surface temperature T = 250 K and radius R = 1km. 

Assuming a conical ejection fan of width 120° and a power-law index a = 3.9, we get 
good fits for certain fan directions by using C v = 36cm 1 / 2 ms~ 1 . An example of a model 
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spectrum giving a good fit is shown in Fig. 3. This example corresponds to an Earth- 
comet-fan direction of (165°, 0°). The coordinate pair represents cometocentric azimuth and 
elevation relative to the comet orbit plane, with azimuth given relative to the Sun direction, 
so that an antisolar tail would be in the (180°, 0°) direction. The most sunward fan direction 
giving a good fit is (—45°, —15°). A precisely sunward (0°, 0°) fan does not fit the data 
because of its positive Doppler offset, as shown in Fig. 3. The relationships between the Sun, 
the comet fan direction and the viewing geometry are illustrated in Fig. 4, which shows the 
geometry for the three fan directions discussed. 

While we have discussed only two examples, there is a locus of fan directions that give 
a good fit corresponding to Sun-comet-fan angles varying from 50 to 165°, depending on 
the plane-of-sky position angle. The radar observes the radial component of velocity of 
the grains with respect to the Earth. This radial component depends on the plane-of-sky 
direction because of the Keplerian motion of the grain particles. However, good fits are 
obtained for all directions with the same Earth-comet-jet angle, 100° from the Earth-comet 
direction (as shown in Fig. 4). 

The velocity scale factor C v = 36cm 1//2 ms" 1 corresponds to V t = 36 ms -1 for an 
a = 1 cm grain (neglecting gravity). Changing the size distribution to an a = 3.5 power 
law reduces the spectrum width to half that for the steeper power law because it weights 
the larger, slower grains more heavily. For this shallower power law, simply increasing C v 
is ineffectual at widening the spectrum, as a m increases at the same time. Hence, a simple 
gas-drag model suggests a steep (a ~ 4) size distribution near the cm size scale in order to 
match the observations. Although the overall dust size spectral index a averaged about 3.5 
or 3.6 for Hyakutake and Halley, it can be closer to 4 at the largest grain sizes or when the 
grains are undergoing fragmentation (Fulle et al., 1995, Fulle et al., 1997). 

The 36 cm 1 / 2 ms -1 scaling factor (for a = 3.9) is similar to that estimated for the coma 
echo from Comet Hyakutake and significantly larger than that estimated for IRAS-Araki- 
Alcock and Halley, as shown in Table 2. This high C v implies that the gas drag effect for 
C/2001 A2 is strong and capable of lifting off large grains of the size needed to account for 
the observed depolarization. For example, if one assumes R = 1 km and p n = 1 gem -3 , then 
using this C v value in Eq. 4 gives a m = 10 m, which is more than enough to accommodate 
a significant population of large, depolarizing grains. One can see why IRAS-Araki-Alcock 
would have a much smaller a m comparable with the Rayleigh transition size X/2tt, since the 
smaller C v and larger radius for that comet (~ 10 km) would make the C 2 /R 2 factor in Eq. 4 
about a factor of 100 smaller than for C/2001 A2. Radio observations of OH molecules give 
line widths suggesting that C/2001 A2 has a comparable gas outflow velocity (about 600 
ms -1 ) to Halley and Hyakutake (Lovell et al., 2005, Crovisier et al., 2002). Thus, nucleus 
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size may be the most important contributor to large-grain ejection. 

If one inserts C v = 36 cm 1 / 2 ms _1 into Eq. 3 and assumes C D = 2, p g = 0.5 gem -3 , and 
a grain temperature of 250K (implying a gas velocity of 270 ms" 1 ), then solving for Z gives 
a value of 7 x 10 -4 g cm~ 2 s _1 for the gas flux, which is similar to that obtained for Hyakutake 
(Harmon et al., 1997) and an order of magnitude larger than the nominal sublimation rate for 
clean water ice at 1 AU. Such high apparent gas fluxes might be possible if the grain ejection 
were explosive or jet-like, also consistent with the observed OH line width velocity of 600 
ms^ 1 Lovell et al. (2005). A high Z would also give a more reasonable surface active fraction 
than the extremely high values obtained when one combines the nucleus size upper limit and 
measured gas production rates for this comet with the canonical Z value for sublimating ice. 
On the other hand, the observed velocities could be achieved with a smaller Z if the grains 
were extremely fluffy (with a higher drag coefficient and lower grain density). Also, the gas 
flux at the nucleus would be irrelevant if rocket forces from their own outgassing accelerated 
the grains (as discussed later). 



4.2. Mass Loss 

Our basic model of stable and continuous grain ejection can also be used to estimate 
the rate of mass loss M from the central (nuclear) source. This involves calculation of the 
grain production rate required to replenish the population of grains traversing and ultimately 
exiting the radar beam. One approach is to calculate M as a function of the gravity-limited 
size a m assuming C v to be an implicit function of a m via Eq. 3, as was done for IRAS- 
Araki-Alcock (Harmon et al., 1989). Alternatively, one can estimate C v from the spectrum 
width and then estimate M as a function of some non-gravitational size cutoff a max (< a m ), 
as done for Hyakutake (Harmon et al., 1997, 1999). For C/2001 A2 we take the latter 
approach, given that our estimated gravitational cutoff a m from Eq. 3 is much larger than 
A/2-7T and therefore uninteresting as a constraint on M. In this case, assuming isotropic grain 
ejection, the mass-loss rate required to give an observed radar cross section a is given by 



M = a 

where 



8C v p g 
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(4 - a)nhl 



I = L d-a/a m) ^ da ' < 6 > 
Here h is the half-width of the radar beam at the comet and Qb is the grain backscatter 
efficiency computed from Mie theory assuming spherical grains. In Fig. 5 we show the 
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calculated M curve for a = 3.9 and a Q = 1/im. We also plot a curve for a Q — 1 mm in 
order to show the M contribution from the largest grains. The steep slope at low a max 
reflects the strong Rayleigh dependence for Qb- However, since we know that a max > 2 cm 
from the depolarization, then M must actually lie somewhere on the flat part of the curve. 
This means that the large-grain M was ~ l-3xl0 6 gs _1 . For comparison, mass-loss rates in 
(1-10 //m) dust from this comet were estimated from optical spectra to be ~ 3 x 10 5 gs _1 
between late June and mid- July, 2001 (Schleicher and Greer, 2001; Rosenbush et al., 2002), 
while gas production rates during this same period varied between l-2xl0 6 gs _1 during 
quiescent conditions (Schleicher and Greer, 2001; Lecacheux et al., 2001) to 6-9xl0 6 gs _1 
during outburst (Feldman et al., 2002). Thus, we find that the mass in large grains is about 
10 times the mass in 1-10 /im dust, and comparable with the quiescent gas production. 
Table 2 shows that all of the comets with observable grain coma show M ~ 10 6 gs~ 1 . It is 
important to note that our estimated mass-loss rates assume that the grains remain intact 
from the time they leave the nucleus to the time they exit the radar beam. If the grains 
were, in fact, disintegrating or evaporating during this time, then the beam would not be 
filled and our M curves could underestimate the true mass loss and large grain production. 
Conversely, if there is a continuous resupply of grains from the breakup of larger boulders or 
nucleus fragments, then the calculated M is an overestimate. The beam diameter was about 
16000 km for these observations (see Table 1). For comparison, the IRAS-Araki-Alcock radar 
echo was apparently contained within 2000 km (Campbell et al., 1983, Harmon et al., 1989), 
but as that comet showed no depolarized echo, it presumably had smaller grains that could 
evaporate more quickly. 



5. Fragmentation and Grain Production 

The comet was observed to fragment several times before and after we observed it 
(Sekanina et al., 2002). Although our assumption so far has been that the large grains were 
gas-drag ejecta from the surface of nucleus B, the fact that C/2001 A2 was a fragmenting 
object raises other possibilities. An interesting alternative is that the large grains were 
produced by, or were secondary products of, one or more of the observed episodes of nucleus 
fragmentation or splitting. Six fragments (A, G, C, D, E, F) were observed to separate from 
the main object (B) between late March and mid- June, 2001, and three visual- magnitude 
outbursts (presumed to be associated with some of the fragmentation events) were seen on 
March 28, May 10, and June 5 (Sekanina et al., 2002). A fourth outburst with no identifiable 
fragment association was observed on July 11, two days after the last radar observation. 
The possibility that the radar-reflecting grains were related in some way to the last (D-E-F) 
fragmentations of June 7-11 is consistent with the antisolar motion of these fragments. In 
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fact, the observed 210-240° position angles observed for these separating fragments (relative 
to the main nucleus B) between mid and late June (Sekanina et al., 2002) also fit the radar- 
observed velocity, projected from its most antisolar possible direction, 165° from the Sun. 
Grains accelerated to the observed velocities during these fragmentations would have cleared 
out of the antenna beam by the time of the radar observations: Extrapolating the model 
of Sekanina et al. (2002) to our observation dates, fragment D should have left the beam, 
but fragments E and F may have remained within it. Fragments D-E-F had disappeared 
more than a week before these observations, but there may have been remnants. Moreover, 
grains and boulders produced in secondary breakups of fragments might tend to remain 
moving at the same slow (~ms _1 ) velocities as their parents. However, if grains produced in 
secondary fragmentation contained volatiles, then the rocket force produced by outgassing 
from their sunward sides could easily accelerate the grains down the antisolar tail to the 
observed velocities (Harris et al., 1997; Hughes, 2000; Bockelee-Morvan et al., 2001). The 
grain velocity from the rocket force after a time t is given by V = 3V g Zt/iiap g . Hence, 
grains can accelerate to the velocities shown in Fig. 2 in as little as an hour, depending 
on the grain volatile content, though it is not clear whether the acceleration is directional 
enough to produce the observed velocities without simply pinwheeling the grains. Even 
this pinwheeling could affect the resulting spectrum, by increasing the apparent velocity 
dispersion of the grains: a 10 cm radius grain rotating at 10 revolutions per second would 
give an echo bandwidth of up to 100 Hz. Pinwheeling would not explain the asymmetry 
in the spectrum, however. It is possible that a train of icy debris left over from the D-E- 
F fragmentation event, and undergoing continuous secondary fragmentation, could account 
for the radar-reflecting grain population observed long after the major fragmentation event 
itself. The gas produced in such secondary fragmentation could also account for the coma 
wings or "arclets" that were observed symmetric about the tail axis of this comet. Jehin et 
al. (2002) reported arclets from observations on May 16 and July 13. In addition, Woodney 
et al. (2002) observed CN arcs on June 29-30 that were symmetric about a 250° position 
angle axis. Such arcs, which were also observed with comets Hyakutake and C/1999 S4 
(LINEAR), may be produced when gas from a tailward debris train interacts with gas from 
the main nucleus (Harris et al., 1997; Rodionov et al., 1998; Boehnhardt, 2002). The possible 
analogy with Hyakutake is particularly relevant here, as Harmon et al. (1997) suggested that 
the radar coma echo from that comet may have come from the condensation of large grains 
invoked by Harris et al. (1997) to explain the gas arcs and other features. 

It is possible that the radar coma grains were not products of the D-E-F fragmentation 
but rather of some lesser fragmentation event that occurred sometime between the D-E-F 
event and July 7 and which did not produce any large visible fragments. For example, the 
CN arcs during June 29-30 (Woodney et al., 2002) could be an indicator of a late- June 
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fragmentation event, although they may also have been a delayed response to the D-E-F 
fragmentations. Another interesting possibility is that the radar grains were produced in 
some unobserved fragmentation event just preceding and triggering the July 11 outburst 
(labeled "outburst IV" by Sekanina et al., 2002). This seems plausible given that outbursts 
I, II, and III (which Sekanina et al. attribute to secondary disintegration of fragmentation 
debris into fine dust) all peaked a few days after their associated fragmentation events. 
Hence, fresh fragmentation debris may have been present during the radar observations and 
before the outburst itself. The arclet and narrow tailward dust spike observed by Jehin et 
al. (2002) on July 13 might have been gas and dust products of the fragmentation of the 
same large grains responsible for the radar coma echoes. Since it would have taken at least 
four days for a 1-cm grain to clear out of the radar beam (at the maximum observed velocity 
of 50 ms _1 ), some of the radar coma grains associated with an outburst-IV fragmentation 
could also have come directly from the nucleus rather than from secondary fragmentation. 
The fact that the days up to and including outburst IV were likely to have been a time of 
chaotic activity suggests that one might have expected to see some day-to-day changes in 
the radar echoes. This may, in fact, explain some of the echo shape variations that were 
seen. Finally, it is worth noting that the observation of the Hyakutake coma echo (Harmon 
et al., 1997) was made only four days after a major outburst from that comet (Schleicher and 
Osip, 2002). Hence, it is possible that the coma echoes for both C/2001 A2 and Hyakutake 
reflected enhanced grain injection associated with outbursts, although the time sequence for 
the radar observations and the outbursts was reversed between the two comets. 



6. Comparison with Other Comets 

Table 3 shows the model results from Section 3, along with total water production rates 
and velocities estimated using other observations. The required gas flux per unit surface area 
Z and particle acceleration parameter C v are similar to that observed in Hyakutake, and 
quite different from that of Halley and IRAS-Araki-Alcock. For this purpose, C/2001 A2 
and Hyakutake have high grain velocities (HGV) and Halley and IRAS-Araki-Alcock have 
low grain velocities (LGV). The model parameter C v is about a factor of 6 higher and the 
parameter Z about a factor of 100 higher in the HGV pair. We can also compare the total 
gas production from other observations. 

The observed gas velocities are all within a factor of 2 of 1 kms -1 , and do not correlate 
with the C v parameter. Unfortunately, the estimate of gas velocity for IRAS-Araki-Alcock 
is very uncertain. The velocities should be considered further with data on more objects. 

The model uses estimates of the gas flux per unit surface area Z. We compare these to 
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the observed gas production rate Q from passive observations of the OH radical at 1667 MHz. 
We use the estimated surface area A of the comets and compute Q/A for comparison. We 
expect Q/A to be smaller than Z, because the comet is not active over its entire surface. 
We can then further compute an active fraction / = {Q/A)/Z. 

Table 3 shows that the HGV objects have areal gas fluxes Q/A that are larger than the 
LGV objects. The difference is perhaps not conclusive for C/2001 A2, but recall that we 
only have an upper limit on the size of its nucleus. Since the area goes as R 2 , the areal flux 
may easily be substantially larger. 

The correlation between gas production (per unit area) and large-grain velocity seems 
significant. They may both be related to the fact that both C/2001 A2 and Hyakutake 
suffered breakup events (Sekanina et al., 2002, Schleicher and Woodney, 2003), which likely 
produced a lot of gas and dust. C/2001 A2 and Hyakutake were also both much smaller 
than Halley and IRAS-Araki-Alcock, which may have influenced dust production as well. 

7. Conclusions 

The radar observations of C/2001 A2 (LINEAR) place an upper limit of only 3 km for 
diameter of the nucleus and yet yield a grain-coma echo with a cross section second only 
to that from Comet Halley. These results, combined with other (non-radar) observations, 
establish C/2001 A2 as a small but active comet similar to Hyakutake. Like Hyakutake, this 
comet showed characteristic grain velocities of several tens of meters per second, substantially 
higher than for the grain-coma echoes from IRAS-Araki-Alcock and Halley. 

Perhaps the most important new result from these observations is the detection of 
significant depolarization in the coma echo. The earlier radar observations of comets Halley 
and Hyakutake showed hints of depolarization, and mass-loss arguments had always implied 
the dominance of cm-size or larger grains in the radar coma echoes. Nevertheless, this solid 
detection of substantial echo depolarization provides the first unambiguous radar evidence for 
coma grains larger than the Rayleigh transition size X/2tt = 2 cm. The marked contrast with 
the IRAS-Araki-Alcock coma echo, which showed only 1.4% depolarization, is explainable 
with the simple gas-drag theory. The higher gas flux (inferred from the high grain velocities) 
and smaller nucleus size for C/2001 A2 would be expected to give a gravitational grain 
size limit much larger than the Rayleigh transition size. The fact that C/2001 A2 was a 
fragmenting and outbursting object suggests that a simple gas-drag model with a nucleus- 
centered source may not tell the whole story for this comet. The observation of antisolar 
fragmentation, gas arcs, and dust trains suggest that secondary grain fragmentation and 
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tailward self-acceleration could have been important processes. These processes are at least 
consistent with the negative Doppler offset of the echo and could also explain the high grain 
velocities. We note that these same comments could apply just as well to Comet Hyakutake, 
another fragmenting and outbursting comet showing similar radar characteristics. Thus, 
while it is not clear that the simple gas-drag model adequately describes C/2001 A2, it 
provides a framework in which we can determine that the radar data require a reasonably 
compact fan of ejected grains in a direction no closer than 45° to sunward, and also a steep 
(radius exponent a ~ 3.9) size distribution of particles at cm-scales. Significant mass loss 
for this comet occurs in the form of cm-size grains, suggesting that large grains may also be 
present in other small, active comets. Future comet radar observations and comparisons with 
optical and other measurements should shed more light on large-grain production processes 
in active comets. 
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Table 1. Observing circumstances and radar measurements. 



UTC Mid-Time 


RTT 


Runs 


A 


R 


Beam 


Angle 


Integration 


(Toe 




Total a 


He 




00 




(AU) 


(AU) 


(km) 




(s) 


(km 2 ) 


(km 2 ) 


(km 2 ) 




2001 July 7 08:53:44 


262 


13 


0.262 


1.132 


16155 


58° 


3210 


4.10 


1.20 


5.30 


0.27 


2001 July 8 09:29:53 


268 


4 


0.268 


1.145 


16497 


55° 


920 


3.36 


0.82 


4.18 


0.24 


2001 July 9 08:12:14 


274 


8 


0.274 


1.158 


16878 


53° 


1920 


4.98 


1.47 


6.45 


0.30 


Sum 














6050 


4.43 ± 1.33 


1.13 ±0.34 


5.56 ± 1.67 


0.28 ±0.03 



Note. — Echo power results arc given for each day of observations and the variance-weighted sum of all days. RTT is the round-trip light time to 
the object. "Runs" is the number of transmit-receive cycles. Each cycle is 2 RTTs long. A is the Earth-comet distance. R is the comet-Sun distance. 
"Beam" is the diameter to the l/\/2-power points of the telescope beam at the distance of the comet, or half-power considering the two-way path. 
"Angle" is the Sun-comct-Earth angle. <j oc and cr BC are the OC and SC radar cross-sections. /j, c is the circular polarization ratio cr BC /<j oc . The random 
errors (due to thermal noise in the telescope and detector) for the cross-sections are a few percent. However, we estimate the systematic calibration 
uncertainties to be 30%, as reported for the summed results. The uncertainties in the individual values arc discussed in the text. The day-to-day 
variations in cr oc are > 20%, and probably reflect real changes in the comet. 



Table 2. Parameters for all comets for which grain coma was detectable by radar. 



Comet 


Date 


A 


coma 


He 


2R 


a 




Cv 


Z 


M 






(AU) 


(km 2 ) 




(km) 






(cm 1 / 2 ms -1 ) 


(g cm -2 s -1 ) 


(gs- 1 ) 


C/2001 A2 


2001-07-08 


0.27 


4.43 


0.28 ± 0.03 


2 


3.9 


10 m 


36 


70 x 10~ 5 


1-3 xlO 6 


Halley a 


1985-11-29 


0.63 


32 


0.52 ± 0.26 


10 


3.5 


> 2 cm 


5 C 


0.8 x 10- 5b 


0.5-2 x 10 6 


IRAS-Araki-Alcock c 


1983-05-11 


0.03 


0.80 


0.014 ±0.003 


8.8 


3.5 


3 cm 


8 


1.2 x 10~ 5 


0.3-1 x 10 6 


Hyakutakc d 


1996-03-24 


0.11 


1.33 


0.31 ± 0.12 


2.5 


3.5 


> 1 cm 


40 


40 x 10~ 5 


1 x 10 6 


C/2002 06° 


2002-08-09 




1.1 


0.32 ± 0.08 















Note. — The dates are approximate mid-times. A is the Earth-comet distance at the approximate mid-times. The coma cross sections all 
have systematic calibration uncertainties of 30%. 2R is the comet diameter assumed in the modeling, —a is the power-law exponent for the 
grain size distribution. C v is the velocity scale factor, either from a model fit or from Eqs. 3 and 4. Z is the gas flux at the surface required to 
lift particles of size a max off the surface. M is the estimated mass in grains required to match the observed radar cross section using the model 
size distribution. Campbell ct al. (1989) suggested a m for Halley may actually be larger than 2 cm, giving a proportionately larger Z. The same 
argument applies to Hyakutake. 

a Campbcll ct al. (1989). The low SNR prevented a reliable SC detection. 

b C v and Z for Halley are computed from a max and 2R using Eq. 3 and 4, T = 200 K, p n = 1, p g = 0.5, and Cd = 2. 
c Harmon ct al. (1999) 
d Harmon ct al. (1997) 

° Harmon ct al. (2004). No modeling has been done for C/2002 06. Its spectrum is qualitatively similar to the spectrum of C/2001 A2, with a 
similar bandwidth, but the SNR is lower. 
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Table 3. Model and measured gas parameters. 



Comet 






Z 


Q 


Q/A 


f 


^h 2 o 




(AU) 


(cm 1 / 2 ms _1 ) 


lO^cm-Ss" 1 ) 


10 28 (s- 1 ) 


10 18 (cm- 2 s- 1 ) 




kms -1 


C/2001 A2 


1.14 


36 


23.4 


0.9 


0.03 a 


0.005 a 


0.6 b 


Hallcy 


1.50 


5 


0.268 


4-5 c 


0.016 


0.060 


1.2 


IRAS-Araki-Alcock 


1.004 


8 


0.402 


2.5 


0.01 


0.025 


1-2. 5 d 


Hyakutakc 


1.08 


40 


13.4 


30 c 


1.52 


0.113 


1* 



a Since we only have an upper limit to the diameter, this is a lower limit to the areal flux and active fraction. If the diameter 
is really 1 km, then this flux is 0.30 and the active fraction is 0.05. 

b Lovcll ct al. (2005). 

c Gerard ct al. (1987). 

d Irvine et al. (1984). The reported detections were very low SNR. 
e Gerard et al. (1998). The observations were from 1996 March 24. 
f Lovcll (1999). 

Note. — The gas flux Z is from Table 3, converted to molecules per second of H2O. The total gas production rates Q arc 
measured values from the literature as indicated. Q/A is converted to an areal flux using the size estimate in Table 2, assuming 
all gas is from the nucleus. The active fraction / is estimated as (Q/A)/Z. 



-24- 




Fig. I. — Doppler spectra (OC and SC polarizations) for Comet C/1983 HI (IRAS-Araki- 
Alcock) show the narrow nucleus echo and broad grain-coma echo. The SC/OC polarization 
ratio of the coma is estimated to be /i c = 1.4%, after correcting for instrumental crosstalk. 
The spectrum has been truncated so that only the bottom 2% of the nucleus echo show. 
These data are from Arecibo S-band (A = 13 cm) radar observations on 1983 May 11 (Harmon 
et al., 1989). 
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Fig. 2. — Echo power spectrum of C/2001 A2 on 2001 July 7-July 9, and the variance- 
weighted sum of all three days. The top line shows the OC and the bottom line the SC 
polarization. The SC spectra have been vertically offset by -0.007 km 2 /Hz for clarity. The 
spectra were taken at a resolution of 10 Hz, and then smoothed to 40 Hz. 
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Fig. 3. — Radar Doppler spectrum with an overplotted model spectrum (heavy solid curve) 
for a 120° grain ejection fan centered on the near-antisolar direction (165°, 0°) and assuming 
C v = 36cm 1 / 2 ms" 1 and a = 3.9. A model spectrum for a sunward (0°, 0°) fan using these 
same parameters (short dashes) gives a much poorer fit. 
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(a) 




Fig. 4.— Orbital geometry for C/2001 A2 (LINEAR) during June-July, 2001, as seen (a) 
looking down on the ecliptic plane and (b) in the ecliptic plane. The Earth-comet positions 
are shown for June 15 and at the time of the radar observation on July 7. Also shown are 
July 7 cometocentric vectors showing the sunward direction (0°, 0°) (a), the near-antisolar 
direction (165°, 0°) (b), and the (—45°, —15°) direction (c). Both b and c are consistent 
with the radar data, because they have the same Earth-comet-fan angle of about 100°, but 
a is not. 
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Fig. 5. — Mass-loss rate M (solid curve) computed from Eq. 4 assuming C v = 36 cm 1 / 2 
ms -1 , p g = 0.5 gem -3 , a = 3.9, and a — 1 fim. Also shown (dashed) is the curve computed 
assuming these same parameters but with a Q — 1 mm. 



